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Matter effects on neutrino oscillations in both, a supernova and the 
Earth, change the observed supernova neutrino spectra. We calculate the 
expected number of supernova neutrino interactions for ICARUS, SK and 
SNO detectors as a function of the distance which they traveled in the 
Earth. Calculations are performed for supernova type II at 10 kpc from 
the Earth, using standard supernova neutrino fluxes described by thermal 
Fermi-Dirac distributions and the PREM I Earth matter density profile. 

PACS numbers: 13.15.+g, 14.60.Pq, 97.60.Bw 

1. Introduction 

Iron is the most strongly bound of all elements which means that fusion 
and fission reactions result in the absorption, rather than the production, of 
energy. The last stage of a massive star life (mass > IOMq) is the collapse 
of its Fe core (the whole core material has already been transformed, via 
the chain of nuclear reactions, into Fe). This happens when the mass of 
the Fe core exceeds the Chandrasekhar limit (-^ 1.45Mq if one takes equal 
numbers of neutrons and protons, T) — a supernova (SN) type II birth is 
a fact. It starts to explode. All flavors of neutrinos are radiated away in the 
form of two bursts of the duration of milliseconds and seconds. The emitted 
neutrinos carry almost all (~99%) of the SN binding energy (~10^^ erg). 

The only observed, up to now, burst of SN neutrinos came from the 
SN1987A which had exploded in the Large Magellanic Cloud, at the dis- 
tance of about 52 kpc away from the Earth. Due to this distinct distance, 
the reconstruction of only 19 events of the neutrino interactions by the 
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Kamiokande [2] and 1MB |3] water Cerenkov detectors had been possible. 
But, it was enough to confirm the main features of the models of SN explo- 
sion. 

Neutrinos, on their way from the production point inside the SN high 
dense Fe core to the terrestrial detector, interact with matter. Non-zero 
neutrino masses, together with flavor mixing enhanced by matter effects, 
result in considerable differences in neutrino fluxes between the production 
and the detection points. Systematic studies of the number of SN neutrino 
interactions in three detectors, ICARUS 4_ , SK 5 and SNO ,6 , as a func- 
tion of the distance passed by neutrinos in the Earth is the aim of this 
work. 



A supernova is a source of fluxes of neutrinos and antineutrinos of all 
three flavors (e, r). The neutrino/antineutrino energy spectra at the 
production point {i.e. inside the SN) of various flavors can be described by 
the thermal Fermi-Dirac distributions (all chemical potentials set to zeros. 



where a = v^, f/^, z^r, z^ej Vr- Here, E represents energy of the neutrinos. 
La is the total energy released in various flavors of neutrinos [La Eb/6, 
with i?B — 3 X 10^^ ergs — binding energy emitted in the core collapse of 
the star), is a normalizing constant given by F3 = 77r''/120, Ta is the 
temperature of the I'a gas in the neutrino sphere. We assume the following 
hierarchy of temperatures: T^^ = 3.5 MeV, Tp^ = 5 MeV, T^^^p^ = 8 MeV, 
where Vx and mean f^, z^,- and P^, i/r, respectively. 

Neutrinos produced deep inside supernova, before escaping from the star, 
interact with the matter of its mantle. Thus, transitions of neutrino species 
can occur. These oscillations of neutrinos in supernova are considered here 
according to |S] . The conversions take mainly place in two resonance layers 
in the outer regions of a supernova mantle. There is a high density reso- 
nance region (H resonance layer) — Ph ~ 10^-10'^ g/cm'^, and low density 
resonance region (L resonance layer) — /Ol ~ 10-30 g/cm^. Transitions 
in the H resonance layer are governed by atmospheric neutrino oscillation 
parameters: Am^^ and Ois, whereas transitions in the L resonance region 
are governed by solar neutrino oscillation parameters: Am2i and 6*12. The 
transition probabilities (between two neutrino mass eigenstates) in the res- 
onance layers are called the flip probabilities — Pu, Pl for neutrinos and 
Ph) Pl for antineutrinos. 



2. Neutrino fluxes from a supernova 



see 13): 
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The final total neutrino/antineutrino mass eigenstates fluxes on a super- 
nova surface {Fi,Fi, i = 1,2, 3) are given by the following sets of equations 
{F^, F^ are the initial total production fluxes defined above, respectively, e 
— i^e, e — i^e, X — or Vr, X — i/^ or Ur). 
For Direct (normal) mass hierarchy (mi < m2 <C m^): 

Fi = PuPlF!^ + (1 - PrPl)fS , 

F2 = {PH-PnPL)F^ + {l-PH + PHPL)F^, 

F3 = {1 - Pu)F^ + PnFl^ , 

F-, = (l-PO^o + PLi^o, 

F2 = PlF^ + (1 - Pl)F^ , 
F, = F° 

For Inverted mass hierarchy (711,3 1^1 < "^-2): 

PlF0 + (1-Pl)F°, 
(1-Pl)FP + PlP0, 

^ X 1 

(Ph - PhPl)P° + (1 - Ph + PhPl)P| , 
PhPlP" + (1 - PhPl)P° , 
(l-PH)P| + PHP^ 

A complete discussion of how to calculate flip probabilities can be found 
in '8'. For the purpose of this paper we will only state here that, in case one 
considers the LMA (Large Mixing Angle) neutrino oscillations parameters: 
Pl = Pl = and, in the so called Large O13 case (sin^ O13 > 3 x 10~^, the 
so called region I in 8 ): Ph = Ph = 0, while in the so called Small 6*13 
case (sin^ 013 < 2 X 10~^ the so called region III in [S]): Ph = Ph = 1 (it 
is interesting to notice here that, in this case, the resulting fluxes are equal 
for both. Direct and Inverted, mass hierarchies). 

The finite spread of the neutrino wave packets, together with the small 
value of their coherence length and the large distance from supernova to the 
Earth, imply that neutrinos arrive to the surface of the Earth as fluxes of 
incoherent mass eigenstates. 



Pl = 

F2 = 

F3 = 

F-i = 

F-2 = 

F-3 = 
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Thus, the above equations describe also (except for a simple geometri- 
cal factor related to the distance between a supernova and the Earth) the 
neutrino/antineutrino mass eigenstates fluxes at the surface of the Earth 
(because there is no matter on the way between a supernova and the Earth, 
there can be no additional transitions between neutrino mass eigenstates). 

3. Earth matter effect 

In order to calculate oscillation probabilities we used the standard de- 
scription of neutrino regeneration effect in the Earth P]. The numerical 
calculations were made with use of the CERN library function DEQBS. We 
applied the realistic Earth matter density profile PREM I |101 [llj . 

The following LMA I (Large Mixing Angle) neutrino oscillation param- 
eters were used (note the two values of the sin^ Ois)- 
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In addition we also considered two cases of mass hierarchies, the Di- 
rect (normal, mi < ?7i2 ^ m^) and the Inverted (7713 <^ nii < 771,2) ones 
(however, because the energies of SN neutrinos are relatively small and the 
Dirac's phase 6 is set to zero, there are almost no differences in results of 
the regeneration in the Earth between these two hierarchies). 

4. Detection of supernova neutrinos 

We consider three detectors: ICARUS (U, SK and SNO 0. The po- 
sitions of these detectors on the Earth are shown in Fig. ^ These detectors 
have the following feature: in most of the time, when a possible neutrino 
signal arrives to the Earth, at least one detector is shielded by the Earth 
(assuming a supernova exploded in the center of our Galaxy), and therefore 
we will see the regeneration effect in the Earth. All possible processes which 
contribute to the total number of neutrino interaction in each detector are 
taken into account. The cross sections for these processes which we use 
come from [12' and '13J . 
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Fig. 1. An example of a possible relative position of ICARUS, SK and SNO detec- 
tors and a supernova which exploded in the center of our Galaxy. 

5. Results and discussion 

The expected total numbers of supernova neutrino interactions A'^sn; 
integrated over neutrino energy in the range 0.1 MeV-100 MeV (the whole 
supernova neutrino energy spectrum), for all possible neutrino processes 
in ICARUS T600 (an "industrial" ICARUS module filled with 600 tons of 
liquid argon), SK (32 ktons of light water) and SNO (1 kton of heavy water, 
1.7 ktons of light water) detectors are calculated (neutrino oscillations in 
supernova and the regeneration effect in the Earth have been taken into 
account). It is assumed that a supernova explosion occurred in the center 
of our Galaxy, that is lOkpc away from the Earth. The results are presented 
in Fig.|21as a function of the distance which neutrinos traveled in the Earth, 
for four possible combinations of the mass hierarchy and the 0i3 value. 

The main neutrino interactions with detector materials which contribute 
to the A^sN are the following (the minimum and maximum contributions of 
a particular process, taken for the four considered cases altogether, into the 
total number of interactions are also given): 
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Fig. 2. The expected number of supernova neutrino interactions in ICARUS T600 
(600 tons of liquid argon), SK (32 ktons of light water) and SNO (1 kton of heavy 
water, 1.7 ktons of light water) detectors as a function of the distance, L, neutri- 
nos traveled in the Earth, for different combinations of mass hierarchy and Ois 
(DL — Direct mass hierarchy and Large 6*13, IL — Inverted mass hierarchy and 
Large Ois, DS — Direct mass hierarchy and Small 6*13, IS — Inverted mass 
hierarchy and Small Ois). For details see the description in the text. 
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It can be seen that, while interactions dominate in the ICARUS de- 
tector, z^e interactions dominate in the SK and SNO detectors. Taking into 
account that, in case of Large 6*13, neutrino oscillations in supernova make 
the fe spectrum harder (hot) for the Direct mass hierarchy and the z/g spec- 
trum harder (hot) for the Inverted mass hierarchy (and that all relevant 
cross sections increase with energy), one gets the corresponding behavior of 
A^SN in Fig. [51 (compare the DL versus the IL curves). 

Two conclusions are straightforward: the distance traveled by neutrinos 
in the Earth has only little influence on the value of N^tss and, in case of 
Small 013, the value of Nq,n does not depend on the mass hierarchy at all. 
Finally, the A^sn from all three detectors should allow us to draw conclusions 
about the value of the 0i3 and, in case the 6*13 is sufficiently large, it should 
also be possible to say which mass hierarchy is in force. 

Last, but not least, it should be noted that, in case of the ICARUS 
detector, for the purpose of this paper we performed calculations for the 
currently existing ICARUS T600 (600 tons of liquid argon) module. The 
final total mass of the ICARUS detector (which will be installed in the 
underground LNGS Laboratory in Gran Sasso/Italy, see pO) will be of the 
order of 3000 tons of liquid argon. That means that the expected total 
number of supernova neutrino interactions A'sn will be five times larger 
than the one presented in this paper. 
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